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1.0  INTRODUCTION 

A  frequently  observed  characteristic  of  Synthetic  Aperture  Radar  (SAR)  imagery 
of  moving  ships  is  an  area  of  increased  backscatter  immediately  behind  the  ship  and 
extending  for  approximately  one  ship  length.  The  region  causing  these  signatures  has 
been  termed  the  near-field  disturbed  water  (NFDW)  region  of  the  ship.  As  part  of 
the  Wakes  program  we  have  characterized  the  SAR  signatures  of  this  region  under  a 
variety  of  ship  and  environmental  conditions  (Malinas,  1992),  and  have  begun  a 
modeling  effort  with  the  goal  of  being  able  to  predict  the  appearance  of  these  NFDW 
returns.  The  model  we  have  developed  is  based  on  a  Fourier  transform  solution  of 
the  linearized  boundary  conditions  with  an  applied  pressure  field.  A  procedure  was 
developed  for  calculating  the  pressure  field  for  a  given  hull  shape  and  ship  speed,  and 
this  pressure  field  was  used  for  calculating  the  near-field  Kelvin  wake.  A  wave 
breaking  criterion  was  then  applied  to  the  vertical  acceleration  computed  from  this 
model  in  order  to  estimate  the  extent  of  the  breaking  regions.  The  fluid  velocities 
within  the  breaking  regions  were  computed,  and  the  radar  cross  section  was  assumed 
to  be  proportional  to  the  breaking  area.  The  radar  cross  section  and  fluid  velocity 
were  used  to  generate  a  simulated  SAR  image  which  can  then  compared  with  an 
actual  SAR  image.  A  description  of  the  data  set  used  in  this  analysis  is  given  in 
Section  2.  In  Section  3  a  model  for  the  NFDW  returns  is  presented  and  the  model 
predictions  are  analyzed  and  compared  with  SAR  measurements  in  Section  4. 

Finally,  in  Section  5  conclusions  are  made  and  a  recommended  course  of  action  is 
presented. 
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2.0  BACKGROUND 

SAR  observations  of  range  travdiing  ships  are  useful  for  investigating  the 
NFDW  signature  because  the  ship  is  displaced  from  the  wake  thus  allowing  an 
analysis  of  the  wake  feature  without  significant  interference  from  the  ship  return. 
These  images  are  also  interesting  because  an  azimuthal  displacement  of  die  Mfak&  itself 
has  been  observed,  indicating  that  relatively  large  Doppler  shifts  or  scat^erer  velocities 
occur  within  the  NFDW.  For  the  work  described  in  this  report,  m  r  rticular  data 
set  of  this  type  was  chosen  and  additional  parametric  variations  were  considered. 

SAR  data  were  collected  by  the  ERIM/NAWC  P-3  during  the  1989  OMR  Ship 
Wake  Experiment  (Lyden,  1989)  under  a  range  of  operating  and  imaging  conditions. 
The  SAR  imagery  for  FFO  Run  7-2  (Figure  1)  was  collected  with  the  ship  traveling 
in  the  range  direction  at  17  knots  and  at  a  relatively  large  incidence  angle  (65 
degrees).  In  this  data  the  ship  is  displaced  in  azimuth  due  to  the  doppler  shift 
associated  with  the  ship’s  velocity.  Bright  returns  arc  located  in  range  near  the  bow 
and  stem  of  the  ship  extending  in  azimuth  from  the  actual  ship  position  toward  the 
displaced  ship  position.  Thus,  these  returns  appear  to  have  dealer  shifts  ranging 
from  nearly  zero  up  to  a  maximum  value  which  is  near  the  ship  q;)eed.  Furthermore, 
the  features  are  slightly  tilted,  indicating  that  the  doppler  shift  appears  to  decrease 
with  distance  aft.  Measurements  of  the  apparent  wake  velocity  are  presented  in 
(Malinas,  1992)  for  this  particular  case  and  other  range  traveling  cases  as  well. 

These  measurements  correspond  to  the  maximum  observed  velocities  of  the  wake 
feature  for  each  range  cell. 

Because  of  the  range  locations  of  these  bright  returns,  it  seems  reasonable  to 
associate  them  with  the  breaking  bow  (or  shoulder)  and  stem  ^  /aves.  Figure  2 
contains  aerial  photos  of  frigates  traveling  at  two  speeds  illustrating  these  breaking 
waves.  The  return  near  the  stem  is  brighter  and  extends  to  larger  Doppler  shifts  than 
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Ship  Speed:  17  knots  Incidence  Angle:  65  deg.  Wind  Speed:  10.2  knots  Rel-  Wind.  Dir:  -84  deg. 
Figure  1:  Four  Channel  Lnoageiy  for  FFG  Run  7-2  1/31  (Range  Traveling  Ship). 


Aerial  Photograph 

FFQ  - 10  (Keel  Dome) 


Figure  2a:  Aerial  Photograph  of  FFG-10  (Kee!  Dome)  at  18  and  25  Knots. 
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Figure  2b:  Aerial  Photograph  of  FFG-10  (Bow  Dome)  at  1 
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the  return  near  the  bow.  We  have  undertaken  a  modeling  effort  in  order  to  try  to 
predict  at  least  some  aspects  of  the  observed  signature.  The  region  near  the  ship  hull 
is  extremely  complex  and  difficult  to  model  adequately,  and  we  realize  that  the  work 
presented  here  represents  a  simplified  and  perhaps  overly  simplistic  treatment  of  the 
problem.  However,  we  hope  that  this  preliminary  study  will  provide  some  basis  or 
motivation  for  additional  investigations  in  this  area. 


3.0  WAKE  MODEL 

The  model  used  in  this  study  is  based  on  a  Fourier  transform  solution  of  the 
linearized  boundary  conditions  with  an  applied  pressure  field 

iK*0^)-F*+ 

in  a  coordinate  system  moving  with  the  ship,  where  the  exponential  growth  term  is 
used  as  a  device  to  satisfy  the  radiation  condition  (Whitham,  1974).  The  velocity 
potential  is  assumed  to  be  of  the  form 

(2) 

where  is  the  ship  speed  (in  the  +x  direction)  and  4  is  the  perturbation  potential 
induced  by  the  ship.  The  linearized  boundary  conditions  are  then 

at  z«0.  The  velocity  potential  must  also  satisfy  Laplace’s  equation,  of  course.  A 
Fourier  transform  solution  of  this  set  of  equations  can  be  written  as 

and 

where  k^^k^+k^.  Substituting  these  into  the  boundary  conditions  yields 


(k,U,*i€f-gk 
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and 


where  Fik^  is  the  Fourier  transform  of  Ax,y). 


(7) 


3. 1  NUMERICAL  KELVIN  WAKE  CALCULATIONS 

During  the  first  phase  of  this  project,  we  used  a  very  crude  approximation  for 
the  pressure  distribution,  namely,  the  hydrostatic  pressure  on  a  rectangular  hull. 

These  results  provided  encouraging  agreement  between  the  modeled  locations  of  the 
breaking  regions  and  the  observations  from  the  aerial  photographs.  During  the 
second  phase  we  mnde  two  improvements.  First,  we  used  a  better  approximation  to 
the  hull  shape  by  digitizing  an  aerial  photograph  of  the  ship  and  applying  a  taper 
along  and  across  the  hull.  Secondly,  we  implemented  a  procedure  for  calculating  the 
applied  pressure  distribution  required  to  produce  a  surface  distplacement  corresponding 
to  the  assumed  hull  shape.  Figure  3  outlines  the  algorithm  implemented  for  this 
modeling  activity.  The  complete  end-to  end  simulation  is  performed  through  three 
programs;  wakel,  wake2,  and  sarwakesim.  This  modeling  procedure  is  as  follows. 

The  surface  elevation  is  calculated  for  a  delta-function  pressure  distribution 
using  the  equation 


*^^dkjk 


(8) 


where  k^*it/Ax,  k^^njiiy,  and  Ax  and  Ay  are  the  grid  spacings  in  the  along- 

track  and  across-track  directions.  The  phase  term  involving  was  added  to 

compensate  for  a  shift  in  the  surface  elevation  pattern  which  occurs  when  the  wake  is 
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computed  on  a  finite  grid  and  with  a  finite  value  of  e  (see  below). 

The  surface  devation  is  specified  at  N  points  on  the  ship  hull  and  the  pressure 
at  these  locations  is  obtained  by  solving  the  set  of  N  equations 

P8m-l 

for  the  pressure  .  This  set  of  linear  equations  is  solved  by  re-writing  Eq.(9) 
as  a  matrix  equation  involving  the  point  source  devation  matrix,  h,  and  inverting  the 
equation  by  multiplying  both  sides  by  A'*,  i.e. 

^  (10) 


where  the  j*  dement  of  hf  and  ^  corre^nds  to  the  location  (x^  ,y^)  on  the  ship 
hull,  and  is  the  surface  devation  at  the  location  (Xy,yp  due  to  a  point  source  at 
the  location  (x^y^).  The  solution  of  Jb'*  is  calculated  using  LU  decomposition  with 
forward  and  back  substitution  (Press  et.  al.,  1986).  Note  that  for  1/^-0,  k  becomes  a 
diagonal  matrix  and  pf  reduces  to  the  hydrostatic  pressure  on  the  hull. 

The  Fourier  transform  of  the  pressure  distribution  obtained  by  this  procedure 
was  calculated  and  the  result  was  used  to  calculate  the  surface  devation  and  tlte 
velocity  potential  function  outside  the  hull  region,  i.e. 


n(v)-V/ . *A 


ai) 


and 


10 


where 


Ocja.^itf-gk 


*^:r**> 


(12) 


The  fluid  vdodty  and  the  apparent  vertical  acceleration  in  an  earth-fixed  coordinate 
system  were  then  calculated  from  the  velocity  potential  as 

and 


(15) 


The  wake  obtained  from  this  procedure  depends  on  the  hull  shape  or  pressure 
distribution,  of  course,  but  also  on  the  value  of  the  parameter  e .  In  the  analytic 
theory,  this  parameter  is  allowed  to  approach  zero,  but  in  the  numerical  procedure  a 
finite  value  must  be  used  in  order  to  avoid  wraparound  in  the  FFT  dontain.  From 
equations  (8),  (11)  and  (13),  the  energy  spectrum  is  equal  to  the  product  of  the  energy 
spectrum  of  the  impulse  response  function  h(,x;y)  and  the  q)ectrum  of  the  pressure 
dii.’^tribution  p(xiy) .  The  spectrum  of  h(x;y)  is  confined  to  a  small  region  surrounding 
the  hourglass-shaped  curve  defined  by 

which  has  a  minimum  wavenumber  (where  ly-0)  given  by  .  Note  that 

Eq.(16)  is  for  a  ship  traveling  in  the  x  direction.  The  spectrum  of  hixjf)  will  rotate 
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as  the  ship  direction  changes.  Figure  4  illustrates  the  spectral  shape  of  a  Kelvin  wave 
and  its  velocity  dq)endence  (again  for  a  ship  traveling  in  the  x  direction). 

The  spectral  energy  density  along  this  curve  is  given  by 

^ 07) 

and  the  3-dB  width  of  the  spectrum  in  the  ib,>direction  is  approximately 

Ml.  4eJk* 

- r — r  •  (18) 


The  total  energy  in  the  spectrum  of  h(x^)  is  then  q^roximately  given  by  the  product 
of  the  peak  spectral  density  and  the  spectral  width,  integrated  in  the  direction. 

This  product  turns  out  to  be  indqiendent  of  the  ship  iqNsed,  provided  that  we  choose: 

(1)  dk^-2AJb^'-4ff/Z.,,  in  order  to  have  two  discrete  wavenumber  samples  within  the 

support  region  of  ; 

(2)  where  in  order  to  keq>  the  same  i:umber  of  transverse 

Kelvin  wavelengths  within  the  computtdional  domain; 

(3)  k^"k^lO  in  order  to  keep  the  spectrum  within  the  FFT  domain,  i.e.  to  have  an 
adequate  sample  spacing  relative  to  the  transverse  Kelvin  wavelength. 

The  spectrum  of  p(,x^)  has  two  components;  one  corresponding  to  the 
hydrostatic  pressure  or  hull  shape,  which  is  assumed  not  to  change  with  ship  speed, 
and  the  other  due  to  the  dynamic  pressure  which  increases  roughly  as  the  square  of 
the  ship  speed.  The  spectrum  of  the  hydrostatic  pressure  is  expected  to  have  an 

oscillatory  structure  (depending  on  the  hull  size)  superimposed  on  a  roughly  to 
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k~*  falloff  (depending  on  the  hull  shape).  For  example,  a  rectangular  hull  with 
dimensions  D,  and  would  produce  a  hydrostatic  pressure  q)ectrum  of  the  form 

“  sincHkJ>J2)sine\k^^)  (19) 

which  falls  off  like  k"^,  while  a  triangular  hull  with  the  same  dimensions  at  the 
waterline  would  produce  a  hydrostatic  pressure  spectrum  of  the  form 


\P,(k^\^  oc  sine*(kJ>JA)sUte\k/>^4) 


(20) 


which  falls  off  like  Since  the  wavenumber  location  of  the  support  region  for 
fiOcJk^  scales  as  glUg,  the  wavenumber  falloff  of  the  hydrostatic  component  implies 

that  the  wake  energy  generated  by  this  component  is  proportional  to  U*  for  the 

rectangular  hull  and  U*  for  the  triangular  hull.  The  wave  drag  Z>^,  which  is  given 
by  the  total  wake  energy  divided  by  the  iength  of  the  computational  domain,  is  then 
proportional  to  and  U,  for  these  two  cases,  and  the  normalized  wave  resistance 
coefficient 


(21) 


is  indq[)endent  of  ship  speed  for  the  rectangular  hull  and  proportional  to  U,  for  the 

triangular  huU.  The  dynamic  pressure,  on  the  other  hand,  is  neariy  a  delta-function  in 
the  spatial  domain  and  thus  its  wavenumber  spectrum  is  nearly  flat.  Since  the 

magnitude  of  this  pressure  spike  is  proportional  to  If/,  its  spectrum  and  the  wave 
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eneisy  produced  by  this  component  are  both  proportional  to  U*,  implying  that  the 
wave  resistance  coefficient  is  independent  of  ship  speed  for  this  component. 

3.2  BREAKING  WAVE  PREDICTIONS 

The  next  stage  in  the  modeling  effort  was  to  use  the  Kelvin  wake  model 
described  above  to  predict  regions  in  the  near  field  region  where  wave  breaking  is 
expected  to  occur.  This  is  a  somewhat  controversial  subject,  and  is  open  to  the 
fundamental  criticism  that  it  is  inconsistait  to  use  linear  wave  theory  to  attempt  to 
predict  the  highly  nonlinear  phenomenon  of  wave  breaking  (this  point  is  addressed  by 
Longuet-Higgins,  1969).  Nonetheless,  encouraging  results  have  been  presented  by 
several  investigators  suggesting  that  the  approach  may  have  merit  despite  these 
concq)tual  difficulties. 

It  is  well  known  that  water  waves  have  a  limited  range  of  amplitudes,  beyond 
which  the  wave  "breaks**  and  the  motion  becomes  chaotic  in  the  region  of  the 
breaking  wave  crest.  Exact  solutions  of  the  wave  equations  exist  for  the  steady-state 
case,  which  show  that  the  wave  crest  becomes  sharper  as  the  wave  approaches  its 
limiting  amplitude.  However,  the  conditions  required  for  the  onset  of  breaking  are 
not  well  understood,  particularly  in  the  general  time-dependent  case  when  a  broad 
spectrum  of  wavelengths  is  present.  Phillips  (1958)  suggested  that  waves  break  when 
the  verticai  acceleration  at  the  crest  exceeds  the  gravitational  acceleration  g,  while 
Longuet-Higgins  (1969)  argued  for  the  value  g/2,  which  is  the  maximum  downward 
acceleration  at  the  crest  of  a  Stokes  wave.  Snyder  and  Kennedy  (1983)  used  this 
ciitedon  to  predict  the  occurrence  of  whitecaps  and  found  reasonable  agreement  with 
limited  field  data. 

An  empirical  wave  breaking  criterion  was  determined  from  a  set  of  laboratory 
measurements  by  CX;hi  and  Tsai  (1983),  and  these  results  were  reinterpreted  by 
Srokosz  (1986)  as  being  consistent  with  a  limiting  dowtiward  acceleration  of  0.4g. 
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However,  it  was  pointed  out  by  Longuet-Higgins  (1983)  that  for  Mte-amplitude 
waves  there  is  a  distinction  between  the  apparent  acceleration,  as  measured  by  a  wave 
probe  at  a  fixed  location,  and  the  real  or  Lagrangian  acceleration  experienced  by  a 
given  parcel  of  water.  In  addition,  he  i^owed  that  the  apparent  vertical  acceleration 
is  theoretically  unlimited  but  the  real  acceleration  for  steady  waves  is  always  less  than 
0.4g.  However,  he  also  argued  that  for  a  superposition  of  waves  with  different 
wavelengths  the  real  downward  acceleration  could  conceivably  attain  a  value  of  g. 
Observations  by  Ewing  et  al.  (1987)  showed  a  maximum  real  downward  acceleration 
of  about  0,38g  and  apparent  downward  accelerations  as  large  as  l,6g,  but  there  was 
no  indication  as  to  whether  these  waves  were  breaking, 

For  our  application,  the  apparent  vertical  acceleration  was  computed  by  noting 
that  in  an  earth-fixed  coordinate  system  the  velocity  potential  function  becomes 

<22) 

so  that  the  vertical  acceleration  in  this  coordinate  system  is  given  by 

This  equation  was  used  to  calculate  the  vertical  acceleration  at  each  point  in  the 
scene,  and  a  threshold  was  applied  to  select  the  pixel  locations  where  wave  breaking 
is  predicted  to  occur  on  the  basis  of  this  criterion.  In  the  following,  unless  otherwise 
noted,  a  threshold  of  g/2  is  used  as  the  breaking  criterion. 

The  effect  of  breaking  is  to  produce  foam  and  small-scale  surface  roughness, 
both  of  which  are  assumed  to  decay  exponentially,  with  different  time  constants.  The 
reflectivity  resulting  from  the  production  and  exponential  decay  of  foam  and/or 
surface  roughness  can  be  modeled  as 
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(24) 

where  i  and  /  are  position  indices  in  the  along-track  and  across-track  directions, 
respectively,  r,  is  the  reflectivity  of  actively  breaking  regions,  is  a  breaking 
index  which  is  equal  to  1  if  ai,<-g/2  and  0  elsewhere,  and 

(25) 

where  p  is  the  temporal  decay  rate  of  the  foam  or  surface  roughness.  This 
expression  must  be  evaluat.id  by  stq)ping  through  the  grid  from  right  to  left,  i.e.  from 
larger  to  smaller  values  of  i . 

The  reflectivity  of  breaking  waves  seems  to  be  on  the  order  of  0.5  for 
electromagnetic  radiation  with  wavelengths  shorter  than  a  few  centimeters  (Lewis  and 
Olin,  1980;  Jessup  et  eU,  1991).  This  has  also  been  confirmed  by  recent  X-band 
laboratory  measurements  at  the  University  of  Michigan  (Walker  et  ah  1993). 
Therefore,  the  same  equation  was  assumed  to  iq^ly,  with  possibly  different  values  of 
P ,  for  both  visible  and  radar  backscatter.  In  the  simulations  to  be  described,  a  value 
of  /3  «0.5  Hz  was  used. 

3.3  SAR  SIMULATIONS 

SAR  images  are  dependent  on  both  the  reflectivity  and  the  motion  of  the  ocean 
surface.  This  dependence  is  described  by  the  SAR  impulse  response  function,  which 
is  proportional  to  the  finite-resolution  Doppler  plectrum  of  the  backscattered  signal, 
i.e.  the  Doppler  spectrum  as  computed  over  a  time  interval  equal  to  the  SAR 
integration  time  (Hasselmann  et  al,  1985;  Lyzenga,  1986).  Thus,  a  collection  of 
scatterers  located  at  the  along-track  position  x  is  mapped  into  a  region  of  the  image 
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centered  at  the  along-track  position  %'  =  and  fiiving  a  width  8jc  = 

where  R  is  the  range,  V  is  the  platform  velocity,  VJix)  in  the  mean  radial  velocity  or 

the  centroid  of  the  D(^ler  spectrum  for  the  scatterers  at  the  location  x ,  and  is 

the  standard  deviation  of  these  velocities,  or  the  width  of  the  Doppler  spectrum. 

In  order  to  simulate  the  SAR  image  of  the  near-field  disturbed  water,  the 
reflectivity  and  the  surface  velocity  was  computed  at  each  grid  cell  as  discussed  in  the 
previous  sections.  The  mean  radial  velocity  was  computed  as 

V,  “  K,coa(^e  +  Kyiln4iein6  ■*-  K,coa6  (26) 

where  6  is  the  incidence  angle  and  is  the  SAR  look  direction,  defined  such  that 
4i*0  if  the  ship  is  moving  toward  the  SAR,  i.e.  thex-axis  points  toward  the  SAR 
ground  track,  and  if  the  y'-axis  points  toward  the  ground  track.  A  nominal 
value  of  ■■O.S  m/s  was  assumed,  and  the  reflectivity  was  distributed  over  the  SAR 

image  using  the  mapping  procedure  described  above.  For  comparison  with  data 
collected  by  the  ERIM/NAWC  P-3  SAR,  the  Doppler  spectrum  was  cut  off  at  a 
frequency  of  0.92m'‘  in  order  to  account  for  the  effect  of  the  presum  filter.  A  more 
detailed  description  of  the  attenuation  effects  of  the  presum  filter  is  presented  in 
(Malinas,  1991). 

3.4  MODEL  PARAMETERS 

Runs  of  the  model  were  performed  to  observe  the  parametric  dependencies  of 
hull  structure  and  ship  speed  and  for  comparison  to  aerial  photos  and  SAR  data.  The 
hull  structures  used,  although  more  complicated  than  the  rectangular  box  shape 
described  earlier,  are  still  very  simple  structures.  Detailed  hull  structures  are 
available  for  several  ship  types  (Kuhn,  1992  ,  Hoekstra,  1991),  however,  the  simple 
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hull  structures  used  in  this  analysis  ate  adequate  for  observing  first  order 
dq)endencies.  The  hull  structures  are  defined  by  three  characteristics;  (1)  the  hull 
shape  of  the  ship  at  the  water  line  from  an  overhead  view  ;  (2)  the  sub-surface  along- 
hull  shape  and  (3)  the  sub-surface  cross-hull  shape.  The  waterline  hull  shape  was 
obtained  by  digitizing  an  aerial  photograph  of  a  firigate  class  ship  (Figure  S).  Three 
types  of  taper  were  considered  for  the  sub-surface  cross-hull;  rectangular,  triangular, 
and  semi-circular  (Figure  6).  These  tapoa  extend  hrom  bow  to  stem  on  the  ship. 

For  example,  the  triangular  hull  has  a  constant  draft  down  the  middle  of  the  ship  from 
bow  to  stem,  regardless  of  the  cross-track  position.  The  along-huU  taper  was  limited 
to  three  possible  cases;  (1)  no  taper  at  all  and  (2)  a  slight  tapering  at  the  bow  and 
stem  of  the  ship  and  (3)  a  taper  at  the  bow  only.  Modifications  to  the  code  for  other 
structures  is  trivial,  however,  these  simple  structures  provide  a  means  for  observing 
the  effect  of  hull  shape  on  the  amount  of  wave  breaking  predicted  by  the  model. 

Note  that  any  hull  structure  may  be  input  into  the  model  with  the  single  limitation  that 
only  one  elevation  may  be  specified  under  the  ship  hull  for  each  location.  Thus, 
bulbous  bows  currently  cannot  be  implemented.  This  limitation  results  from  the 
algorithm  which  calculates  the  pressure  distribution  and  can  be  avoided  by  specifying 
the  pressure  distribution  directly.  In  the  next  section  an  analysis  of  the  results 
produced  by  this  model  is  provided. 
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Frigate  Hull  Shape  -  Top  View 
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CROSS-HULL  SUB-SURFACE  TAPER 
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Figme  6:  Sub-sar&ce  Lx^Ionoited  in  Modd. 
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4.0  RESULTS  AND  DISCUSSION 

In  this  section  an  examination  of  the  model  results  as  a  function  of  speed  and 
hull  shape  is  presented.  First  an  example  run  illustrating  the  model  outputs  will  be 
shown  in  Section  4.1  and  then  comparisons  between  the  trends  predicted  by  theory 
and  those  observed  in  actual  SAR  images  will  be  made  in  Section  4.2, 

4. 1  MODEL  OUl’PUT  EXAMPLE 

Figure  7  shows  the  output  of  the  program  mkel  (refer  to  Figure  3).  The  hull 
shape  in  this  example  has  a  waterline  outline  of  a  frigate  shape  with  a  constant  S 
meter  draft  at  all  hull  locations.  The  upper  left  image  shows  the  surface  elevation  for 
a  point  source  moving  at  9.S  m/s  and  in  the  upper  right  is  the  calculated  pressure 
distribution  for  the  same  speed.  For  the  most  part,  the  pressure  is  flat  across  the  hull. 
There  is  a  notable  increase  in  the  pressure,  however,  at  the  bow  region  which  appears 
to  correspond  to  the  dynamic  pressure  discussed  in  section  3.1.  The  lower  left  image 
is  the  resulting  surface  elevation  for  the  ship  hull  and  the  lower  right  shows  a  scan 
through  the  ship  center  in  the  along  track  direction.  Within  the  defined  region  of  the 
ship,  the  surface  elevation  values  are  equal  to  the  specified  ship  draft.  Modifications 
to  this  hull  shape  wiU  be  discussed  in  Section  4.2.  In  this  example,  the  sample 
spacing  of  the  domain  was  2.9  m  in  both  the  along-track  and  cross-track  directions. 
Recall  from  Section  3.1  that  the  sample  spacing  in  the  along-track  direction  is 
dependent  upon  the  ship  speed,  however  the  cross-track  sample  spacing  is  arbitrary. 

Figure  8  shows  the  different  output  flies  which  are  generated  by  the  programs 
wake2  and  sarwakesim.  This  result  is  for  a  frigate  water  line  shape  but  this  time  with 
a  triangular  sub-surface  shape  in  the  cross-hull  direction.  Again  the  modeled  ship  is 
traveling  at  9.5  m/s.  The  top  left  is  the  surface  elevation  file  which  is  generated  by 
wakel.  The  top  right  image  is  the  radial  velocity  which  is  derived  from  the  velocity 
potential.  The  lower  left  image  is  a  map  illustrating  the  locations  of  the  breaking 
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Figme  8:  Exanq>Ie  of  Program  Wake2  and 


legions  detennined  by  comparison  of  the  vertical  acceleration  at  each  point  with  the 
acceleration  threshold  g/2.  This  map  can  be  compared  to  the  aerial  photographs  in 
Figure  2  where  many  similarities  exist  between  the  structure  of  the  breaking  region 
and  the  white  foam  in  the  photo.  Two  particular  features  stand  out;  (1)  the  bow  (and 
shoulder)  breaking  regions  and  (2)  the  stem  breaking  region  which  fans  out.  The 
bottom  center  image  is  the  reflectivity  which  is  genemted  by  assuming  a  radar  cross 
section  for  a  breaking  cell  and  applying  an  exponential  decay.  Lastly,  the  lower  right 
image  is  the  SAR  simulation  which  is  generated  by  applying  the  SAR  impulse 
response  function  to  the  radar  cross  section.  Observe  that  the  ship  appears  at  the 
bottom  of  the  image,  displaced  due  to  its  radial  motion.  If  the  ship  were  traveling  in 
the  same  direction  as  the  SAR  platform,  much  smaller  Doppler  shifts  would  result  and 
the  wake  signature  would  be  concentrated  nearer  to  the  ship  track.  Note  that  no  SAR 
modeling  of  the  ship  return  itself  is  talcing  place  other  than  shifting  its  position  based 
upon  its  radial  velocity. 

4.2  PARAMETRIC  MODEL  EXAMPLES 

Up  to  this  point,  we  have  provided  an  overview  of  the  model  execution  and 
outputs  but  have  not  yet  compared  the  results  with  tlieory  and  SAR  observations.  In 
this  section  we  provide  a  few  examples  to  illustrate  the  dependence  of  the  breaking 
regions  on  hull  structure  and  ship  speed  and  make  comparisons  to  SAR  data 
observations. 

Figure  9  shows  a  comparison  of  the  modeled  breaking  regions  and  SAR 
intensity  data  for  the  Mgate  example  shown  in  Figure  1.  Results  for  both  X-  and  L- 
bands  are  shown.  The  determination  of  the  breaking  region  is  independent  of 
electromagnetic  considerations,  but  a  frequency  dependence  is  introduced  into 
sarwakesim  by  the  presum  filter  effect  discussed  in  Section  3.3.  Because  the  Doppler 
frequency  shift  due  to  radial  motion  is  inversely  proportional  to  wavelength,  tlie  ship 
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Figure  9:  NFDW  Model  /  SAR  Data  Comparison  For  FFG  Run  7-2. 
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displacement  varies  along  with  the  huli  shape  in  the  examples  shown.  One  would 
expect  an  increase  in  displacement  to  cause  an  increase  in  the  amount  of  breaking 
independently  of  the  hull  shape,  and  the  two  effects  are  combined  in  these  examples. 

We  next  consider  the  effects  of  ship  speed  on  the  amount  and  distribution  of 
wave  breaking.  Figure  11  shows  an  example  of  the  iirigate  case  where  the  triangular 
cross-hull  taper  and  full  taper  along-hull  were  applied.  Because  we  are  making 
comparisons  at  different  speeds,  the  sample  spacing  normalization  described  in 
Section  3.1  has  been  applied.  Thus,  the  slower  the  ship  speed,  the  smaller  the  sample 
spacing  in  the  along-track  direction  and  the  larger  the  number  of  samples  within  the 
ship  hull.  In  all  examples  a  sample  spacing  of  1  meter  in  the  cross-track  direction 
was  used.  This  figure  shows  tliat  as  the  speed  of  the  ship  increases,  the  breaking  area 
increases  roughly  in  proportion  to  the  ship  speed.  This  dependence  is  in  agreement 
with  L-band  SAR  observations  (Malinas,  1991)  which  show  an  increase  in  the  NFDW 
radar  cross  section  of  about  0.3  dBsm  per  knot.  This  comparison  tends  to  support 
our  model,  at  least  under  assumption  that  the  radar  cross  section  is  proportional  to  the 
breaking  area. 

The  last  hull  characteristic  to  be  considered  is  the  hull  outline  at  the  water 
level.  In  all  tlie  previous  examples,  a  frigate  outline  was  used  for  illustration, 
Modifications  to  this  outline  are  considered  in  this  next  example.  Four  different  hull 
outlines  are  implemented  in  this  current  model.  All  have  the  bow  of  a  frigate  and  a 
stem  which  is  selectable.  The  options  for  the  stem  shape  are  (1)  rectangular,  (2) 
triangular,  (3)  ellipticsl,  and  (4)  circular.  Figure  12  illustrates  these  hull  forms  and 
the  sample  spacing  us<xl  for  the  18  knot  case. 

A  comparison  of  the  model  output  for  three  of  these  hulls  is  shown  in  Figure 
13.  The  predicted  dependence  upon  stem  hull  shape  is  very  strong  and  perhaps 
counterintuitive,  with  the  triangular  hull  shape  producing  more  than  10  times  as  much 
breaking  as  the  rectangular  hull.  This  phenomenon  results  from  the  fact  that  the 
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spectrum  of  the  pressure  distribution  for  the  rectangular  stem  is  concentrated  near  the 
A:, -axis,  and  the  spectral  density  in  the  cross-track  direction  falls  off  more  rapidly  for 

the  rectangular  than  tlie  non-rectangular  shapes,  as  shown  in  Figure  14.  The  higher 
spectral  density  at  large  wavenumbers  for  the  non-rectangular  stem  shapes  is  evident 
in  the  short  wavelength  of  the  breaking  waves  shown  in  Figure  13. 

This  figure  also  illustrates  a  general  shortcoming  of  our  model  which  arises 
from  the  neglect  of  energy  dissipation  due  to  wave  breaking.  Because  there  is  no 
dissipation,  waves  which  break  near  the  hull  continue  to  propagate  and  break  at  larger 
distances  from  the  ship  than  is  physically  possible  or  realistic.  That  is,  if  we 
incorporated  an  energy  dissipation  due  to  breaking  waves,  most  of  the  breaking  would 
be  confined  to  a  region  much  closer  to  the  ship  than  is  indicated  in  Figure  13.  This 

would  also  enforce  a  stronger  wavenumber  falloff,  presumably  as  in  the  wake 
spectrum, 
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Figure  14:  Surface  Elevation  Frequency  Response  Comparison  For  Different  Stem 
Hull  Sh£q}es. 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


A  model  has  been  developed  for  prediction  of  the  near-field  disturbed  water 
region  for  a  traveling  ship.  This  model  utilizes  a  Fourier  transform  solution  for  the 
Kelvin  wake  along  with  a  numerical  procedure  for  calculating  the  hull  pressure 
distribution.  The  model  appears  to  predict  the  breaking  regions  more  or  less  correctly 
for  ship  hulls  having  a  rectangular  st^,  but  overpredicts  the  extent  of  the  breaking 
regions  for  non-iectangular  stem  shapes.  This  limitation  is  apparently  due  to  the  fact 
that  we  have  not  included  energy  dissipation  due  to  wave  breaking. 

We  have  used  the  outputs  of  our  NFDW  model  to  simulate  the  SAR  image  of 
this  region  and  found  encouraging  agreement  between  the  appearance  of  the  simulated 
and  actual  SAR  images.  Furthermore^  the  dependence  of  the  total  radar  cross  section 
of  the  NFDW  on  ship  speed  is  in  agreement  with  observations. 

Possible  improvements  to  this  model  would  include  the  incorporation  of  energy 
dissipation  due  to  wave  breaking,  as  mentioned  above,  and  the  development  of  a  more 
detailed  breaking  wave  reflectivity  model  based  on  laboratory  measurements  of  the 
microstructure  of  breaking  waves.  Additional  exercise  of  the  model  with  mote 
realistic  hull  shapes  is  also  recommended. 
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